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Abstract 

Quantum Fourier transform (QFT) is a key function to realize quan- 
tum computers. A QFT followed by measurement was demonstrated on 
a simple circuit based on fiber-optics. The QFT was shown to be robust 
against imperfections in the rotation gate. Error probability was esti- 
mated to be 0.01 per qubit, which corresponded to error-free operation 
on 100 qubits. The error probability can be further reduced by taking the 
majority of the accumulated results. The reduction of error probability 
resulted in a successful QFT demonstration on 1024 qubits. 

1 Introduction 

Quantum computers have been expected to posses more computational ability 
over conventional computers. Shor's factorization algorithmjj, for example, 
has proved the power of quantum computation, which resolves a composite 
number N {L bits long) into the prime factors in the polynomial time of L, 
whereas all the known classical algorithms require exponential resources. The 
factorization problem is also important from a practical point of view, because 
efficient algorithms can break the RSA public-key cryptosystem. 

Although Shor's algorithm and its related algorithms are well established in 
theory, practical implementations of them are still in their infancy. So far, only 
an NMR-based computer has demonstrated the factorization of 15 (= 3 x 5) 
with seven qubits[2]. Until now, it has been difficult to see how one could 
increase the number of qubits to solve problems of practical significance, such as 



breaking the RSA public-key cryptosystem. One of the main obstacles to realize 
large scale quantum computers is errors in gate operations due to decoherence. 
Theory of fault-tolerant quantum computation|3j shows that arbitrary-length 
error-free computation can proceed with only a polynomial overhead in time and 
space, so long as the error probability per fundamental operation is kept below 
a threshold value. Typical predicted values are in 10~^ to 10~^, which depend 
on the product of qubits and gate operations required in the computation, and 
on the model of decoherence |31 These threshold values are still too low for 
the present technology, so that it is desirable to reduce the requirement in the 
scale of quantum computation. 

Recently, Beauregard showed that 2L + 3 qubits are enough to perform 
Shor's algorithm. He used quantum Fourier transform (QFT) for modular expo- 
nentiation, and a QFT followed by measurement (we will refer it to MQFT after 
Measured Quantum Fourier Transform) on the control qubit[Z||Hl. It would be 
remarkable that MQFT is fault-tolerant [Jj, so that the above scheme not only 
reduces the 

number of qubits but also provides robustness against decoherence. 

One of the main issues in implementing quantum information technology is 
selection of quanta to carry the information. Among a number of the implemen- 
tation proposals, photons have played an important role from the beginning of 
the practical proposal on quantum information technology. Photons have ad- 
vantages to implement qubits: SU(2) space easily implemented by polarization, 
very weak coupling to the environment, and existing single photon measurement 
technique. Moreover, we can utilize fruit of extensive research and development 
efforts of the optical communication industry. Commercially available fiber- 
optic devices enable us to construct efficient quantum circuits that consist of 
one-qubit operations (including classically controlled gates) with a reasonable 
cost. Fiber-optics resolves the mode matching problems in conventional optics, 
and provides mechanically stable optical circuits. It would be worth exploring 
feasibility of the quantum information processing based on fiber-optics. In this 
article, we show an implementation of the MQFT on a circuit built up with 
commercially available fiber-optic devices. We have achieved error-free MQFT 
operation on 1024 qubits, despite the imperfections in the actual devices. This 
shows that photonic devices are suitable to process quantum information. 

In the rest of the article, we will report the MQFT experiment in detail. 
We will review the QFT and MQFT used in the phase estimation algorithm 
in section 2. The implementation of MQFT with fiber-optics devices will be 
given in section 3. We will also show successful result on MQFT over 255 qubits 
and 1024 qubits. In section 4, we will consider effects of the imperfections 
in the experimental apparatus, and show the validity of the majority voting 
by repeated measurements. We will also discuss a direction toward making 
quantum computers expanding the present MQFT implementation. 
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2 Theory 



The heart of the factorization and related quantum algorithms lies in the phase 
estimation algorithm "S" that consists of controUed-unitary operations (c-C/'s) 
and MQFT on the control qubits. The phase estimation problem is given as 
follows: An eigenvalue of a unitary transform U defines a phase ip as U \u) = 
exp [27ri if] \u). Our task is to determine the phase expressed in n bits by ip = 
(pi2~^ + • • • + ipn'^~" — 0.(^1 • • • (fin- When the target qubit state is provided 
by the eigenstate |u) of the unitary transform U, the states after controlled-?/^ 
according to the superposition of all the possible control qubit state |j) is given 

by 



2n/2 ^ ' 2"/'^ ' 

The inverse Fourier transform performs the transform 



1 '"-^ 



2"/2 ^ 

J=0 



\^)\u)-.\q,)\u), (2) 



where denotes a state which provide a good estimation of ip when we mea- 
sure the control qubits. In short, the controUed-C/-' provides an unknown phase 
according to the problem. MQFT then determines the phase to find the solu- 
tion. Though a text bookfTU] shows QFT implemented by a quantum circuit 
of Hadamard gates and controlled rotation gates as in Fig. ^a), an alterna- 
tive form of the MQFT circuit was constructed with only one-qubit gates as in 
Fig. nb)[7]. This is based on the fact that c-[/'s commute with measurements 
when the control qubits are measured in the computational basis. We can thus 
replace the controUed-unitary gates with unitary gates controlled classically by 
the results of the measurements [TU| . Since the latter devices act on only the 
target qubit, they are much easier to realize than two-qubit gates. Though 
Griffiths and Niu described their circuit as semiclassical 7 , it is worth to point 
out that the circuit is fully equivalent to the quantum circuit constructed with 
controUed-rotation gates. 

Parker and Plenio 8 found that the phase estimation using MQFT can be 
operated qubit by qubit with only one rotation at a time, as shown in Fig. ^ic). 
This reflects the fact that the RHS of ^ can be rewritten by a product form 
as 

^ (|0> + e2-2"-V (|0) + ,2..2""V . . . ^ ^2.^^ 

The circuit operates on the target qubits, which are in an eigenstate of U , with 
one control qubit at each step. In some applications, however, it is not easy to 
find the eigenstates (if we know an eigenstate, the problem is already solved.) 
Then, instead of using an eigenstate, we can start with a superposition of the 
eigenstates: 



\^^) ^ ^Y.^s\Us) ■ (4) 



s=0 
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We observe \1) — r Ws) in order-finding algorithm [S]. The phase estima- 

tion algorithm proceeds qubit by qubit as follows. In the first step, the c-U^ 
with respect to the first control qubit 2^^^^ (|0) -I- |1)) transforms the initial 
state to an entangled state 

r-l 

(|0) + |1)) \i'o) = (|0) + exp [^z^^J |1)) c, \us) , (5) 

where we use exp [27ri 2"^^iy9] = exp [ni Here (and hereafter,) normalization 
constants are omitted for simplicity. Then a Hadamard gate transforms the 
control qubit to |0) or |1) according to the bit value ipf^ = {0, 1}, and the total 
state is given by 

io) J2 ^s\us) + \i) ^^1"^)- (6) 

The target qubit states collapse to either of the superposition J2s£{>fi= =o} l^s) 
or J2s&{ip= =i}'^s\us) by the measurement on the control qubit. The former 
consists of the eigenstates belonging to the eigenvalues O.xi ■ ■ -Xn^iO, and the 
latter consists of those belonging to the eigenvalues 0.a;i • • -Xn-il, where the 
bit values Xi ■ ■ ■ Xn-i are arbitrary at this step. In the k th step, c-U^ 
provides the state 

(|0) +exp [27ri0.((9j;_fc+i(y5„_fe+2 • • - (fin] |1)) Cs \Us) 

(7) 

The control qubit state is reduced to 2 (|0) -I- exp[27rj 0.</j^_j,_,_j] |1)) by ap- 
plying the classically controlled rotation Rk defined by 

^'^ = (o exp [-27^^$fc]) ' 

where the rotation angle $fc is determined from the bit values obtained in the 
previous measurements, i.e.. 



fe-i ^ 

^ H ^J+Tfn-k+] + l ■ (9) 



Measurement of the control qubit in the computational basis after Hadamard 
transform again yields a deterministic value of (p^-k+i — fn~k+i- Continuing 
the operation steps bit by bit will provide all the bit values of f with probability 
proportional to \cs\^. The target qubit states collapse to one of the binary 
branches according to the measurement result, and thus the number of possible 
target qubit states is reduced at each step. 
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Figure 1: Quantum circuits for measured quantum Fourier transform, (a): A 
quantum circuit with controUed-rotation gates. (b):A quantum circuit with 
classically controlled-rotation gatesjT]. (c):A serial quantum circuit of phase 
estimation H and M stand for Hadamard gate and measurement, respec- 
tively. Phase of the controlled-rotation gates are represented by the numbers 
7r/2, 7r/4, and 7r/8. Rk is defined by © in the text. 
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3 Experiment 



Fig. Ufa) shows an implementation of MQFT. Qubits are represented by the 
single photon polarization states, where the |0) and |1) states are defined to 
be horizontally polarized and vertically polarized, respectively. A sequence of 
single photon pulses enters the quantum circuit. The input photons are generally 
elliptically polarized according to the bit values of the phase ipi, - ■ ■ ,ipn between 
the basis states (|0) and |1)) to simulate the output of the c-U^s as shown in 
©■ 

The key device of the MQFT circuit is the rotation gate, which is imple- 
mented as an interferometer constructed with polarization beam splitters (PBS) 
and a phase modulator (PM). The PBSl divides the states according to the po- 
larizations, and the PM provides a relative phase shift given by © to the |1) 
state. The two polarizations are combined again by the PBS2, and the polar- 
ization state is analyzed by a half wave plate (HWP) and a PBS. The HWP 
provides phase difference between the polarization. When its axis is tilted 22.5 
deg, HWP transforms |0) ^ (|0) + |1)) /%/2 and |1) -> (|0) ~ |1)) /V2, thus acts 
as a Hadamard gate. The accuracy of the phase shift is determined by the 
precision of the electric pulse applied to the PM. In most pulse generators, the 
precision is limited to three digits, i.e., 8-10 bits. This implies the rotation angle 
$fc in © should be truncated at the m th bit (m < 10). 

The interferometer in Fig. EJa) will work in principle, however, is sensitive to 
disturbance in practice. We employed a fiber loop configuration often referred 
to a Sagnac interferometer as shown in Fig. Olb), where orthogonally polarized 
photons propagate in opposite directions through the same fiber. The two basis 
states are subject to the same additional phase fluctuation, so that the present 
MQFT circuit is robust to disturbances. The fiber-loop rotation gate consisted 
of a circulator, a polarization controller (PC), a PBS, a PM, and polarization 
maintaining fiber (PMF). The circulator is a unidirectional device that allows 
light propagation only in the 1^2 and 2 — )■ 3 direction. Before entering the 
PBS, photons were passed through PCI to compensate for the unintentional 
polarization rotation due to the birefringence in the fiber. The photons were 
divided by PBSl. The output ports of the PBSl were both aligned to the 
slow axis of the PMF. The horizontally polarized |0) state propagated in the 
clockwise (CW) direction and the vertically polarized |1) state in the counter 
clockwise (CCW) direction. Note that the PBS converts the polarizations into 
the directions of the propagation, and that the qubits were represented by CW 
and CCW pulses with the same polarization in the fiber-loop. 

The PM was placed at an asymmetric position in the fiber loop to provide 
the relative phase. The CW pulses entered the PM earlier than the CCW pulses 
for 20 ns in the present setting. The arrival time of the CW pulses was syn- 
chronized with the electric pulses to the PM, so that only the photons in the |1) 
state were subject to the phase modulation. To provide accurate phase shifts, 
the electric pulse should be shorter than the difference of the arrival time. The 
present circuit worked with the electric pulses of 10 ns duration. We used a 
high-speed LiNbOs PM designed for 10 Gb/s operation. The CW and CCW 



6 



pulses were combined by PBSl, and went back to the circulator. Qubits were 
again represented by polarizations after PBSl. PC2 compensated the birefrin- 
gence but also provided Hadamard transformation. This was possible because 
PCs creates a fixed but arbitrary phase difference between two orthogonal po- 
larizations, the axis of which can be set also at an arbitrary angle. 

The polarizations of the output photons from the circulator were analyzed 
by PBS2. Photons were detected by a photon detector. The photon detector 
was made of balanced InGaAs avalanche photodiodes cooled to 173 K to reduce 
the dark counts [T^. We obtained the dark count rate at 6-7x10^'' per pulse 
with the detection efficiency of 13 %. The optical loss in the circuit (8.4 dB) 
came mainly from the insertion loss of the PM (6 dB). The whole circuit except 
the photon detector was constructed on a breadboard of 350 mm x 600 mm, as 
shown in Fig. |3| 

The input polarization state expressed by 

_^ (|0) + e^-"-^" |1)) , • • • , ^ (|0) + e^-o-'^--^" |1)) . (10) 

was produced by a motorized polarization controller (Agilent 8 169 A) from the 
pulses of a 1.55-/Ltm laser diode. The pulse duration was about 30 ps. The light 
pulses were then attenuated to 0.7 photons/pulse to ensure that the present 
circuit worked with a single photon. The output bit values were compared with 
the input ones. They were also used to determine the phase modulation to the 
next qubit. 

Figure 01 shows the output photons from the circuit as a function of the 
drive voltage applied to the PM. The input polarization was fixed to the linear 
polarization tilted at 45 deg. Visibility of the interferometer was 98 %, which was 
determined mainly by the extinction ratio of the PBS. We concluded that the 
phase modulation is linear to the drive pulse voltage, because the interference 
signal was well fitted to the cosine curve. The drive voltage for vr phase shift 
was estimated to be 5.80 V from Fig. 0] The phase errors measured in 24151 
rotations for m=5 are shown in Fig. |S1 The values of cos S were distributed 
between [0.98, 1] and [—1, —0.98], with a mean value of (cos (5) = ±0.9936. 

Figureini^a) shows the results of 21 trials of MQFT on 255 qubits. The inset 
shows the result of each trial. 'Successfully transformed bits' in the figure refers 
to the number of bits for which a MQFT operation was done successfully. We 
succeeded making a MQFT of full 255 qubits in two trials. The distribution of 
the successful bits n obeyed a geometric distribution with the error probability 
per qubit p: 

E{n) ^ {1 - pr-^p. (11) 

The average of the successfully transformed bits was 97, which related to the 
error probability by 1/p in the geometric distribution. This reads the value of 
the error probability p = 0.01. More precisely, we estimated the range of the 
error probability as follows. The probability of a successful MQFTs up to the 
fcth qubit n times in N trials is expressed by 

P,= (^^^yE{kr[l-E{k)f-\ (12) 
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Figure 2: Quantum circuit for measured quantum Fourier transform imple- 
mented by fiber-optics. (a):thc principle of the circuit. (b):practical circuit. 
PBS, PM, HWP, and PC stand for polarization beam splitters, phase modula- 
tor, half-wave plate, and polarization controllers, respectively. 
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Figure 3: Quantum circuit for measured quantum Fourier transform constructed 
on a breadboard. 
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Figure 4: Interference fringe of the MQFT circuit. Polarization of the photons 
are set to 45 deg. The results were well fitted by cosine curves. 
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Figure 5: Distribution of the phase error in the phase modulation. The distri- 
bution was measured in 24151 rotations for truncation m—5. 



Suppose the largest (smallest) number obtained in the experiment is kmax 
(kmin)- The maximum value of the error probability per qubit Pmax is de- 
fined so to yield a less than 5 % probability of successful MQFTs for more than 
kmax qubits Nmax timcs. Similarly, the minimum value of the error probability 
per qubit Pmin is defined so to yield a less than 5 % probability of successful 
MQFTs for less than kmin qubits Nmin times. In the present experiment, we 
obtain {kmax^Nmax) = (255,3) and {kmin,Nmin) = (9, 1). The estimated error 
probability per qubit lay in the range 1.2 x 10~^ < P < 2.6 x 10^'^. 

A simple way to reduce error probability is to decide the bit value by majority 
of repeated measurements. We performed a MQFT of 1024 qubit to examine 
the effect on the error reduction by the decision by majority. The rotation 
and measurement were done with the same input polarization states M[— 10) 
times, and the bit value was determined by the majority of the results. We 
set the error probability per qubit dX p — 0.07, intentionally higher than that 
for the previous measurement without accumulation. Figure IHIb) shows the 
results of 30 trials of 1024 qubits. The success probability of full 1024 qubits 
was 80%, and the mean error probability was estimated to be 2.2 x 10^"* per 
qubit. More precisely, we obtained the range of the error probability 1.2 x 10~* 
< p < 4.3 X 10"** with the confidence level of 95 % from the experimental results: 
{kmax, Nmax) = (1024,24) and {kmin, Nmin) = (23,1). The error probability 
satisfies the error threshold for the fault-tolerant operation. 
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Figure 6: Results of MQFT trial, (a) 21 trials of 255 qubits and (b) 30 trials of 
1024 qubits The bit values in (b) were determined by the decision by majority 
of ten measurement outcomes. The inset shows the result of each trial. 'Suc- 
cessfully transformed bits' in the figure refers to the number of bits for which a 
MQFT operation was done successfully. 
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4 Discussion 



4.1 Effects of the imperfection 

In the following, we consider effects of the imperfection of experimental appa- 
ratus. Two types of failure may occur. One is failure in photon detection, the 
other is error in the measurement. If a pulse contains no photons, the operation 
will not affect the target qubits. We can thus continue the calculation by repeat- 
ing the operation step. If a photon is lost in the MQFT circuit, it means that 
the c-U and measurement have been done without knowing the result. There- 
fore, the target states will be in a mixed state corresponding to the two possible 
outcome. Even in this case, we can proceed the calculation by repeating the 
same operation step. The repeated measurement will reduce the target qubit 
state into a pure state by selecting one of the possibilities. 

Errors in the measurement originate from imperfection of the interferometer 
and from errors in the phase modulation. Dark counts are negligibly small in 
the photon detector ^^l- The performance of the interferometer is characterized 
by visibility v, which defines the measurement on the control qubit as 



^^o = a/^|0) (0| + Wi-^|1)(1| 



M, = ^^ |0) (0\ + ^^\l)(l\. (13) 

The phase error, which shifts the rotation angle in JS)) from 27r$fc to 27r$fc+(5, re- 
sults in the control qubit after the rotation gate as 2~^/^ (|0) -I- exp[27ri O-fn-k+i ^ 
The phase error originates from the approximation in the rotation angle $fe and 
from errors in converting the rotation angle into the drive voltage to the PM. 
The latter can be reduced by careful calibration, so that we only have to con- 
sider the effect of the truncation. The truncation at the mth bit results in the 
phase error 

k—l ^ k—1 ^ soo 

= 2^ E <27r5:-^<2vr^-^ = (14) 

j—m j—^n j—^ 

The phase error should not be significant ^Tj, because the contribution from 
the jth bit (j > m) decreases with the factor of 2"^''+-'^^ The worst values of 
cos 5 = ±0.98 obtained in the experiment corresponds to a phase error of 7r/16, 
which agrees quite well with the prediction by H14|l with m—d. The visibility 
and the phase error determine the error probability of the measurement by 

l-wcos^ 
P= 2 ■ 

The estimated error probabilities from H15|l . p — 8.2 x 10^"^ (by using (cos (5) = 
±0.9936) and p = 1.5 x 10^^ (by using cosJmax — ±0.98), agree well with the 
experiment. 
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4.2 Validity of majority voting 



We consider the validity of decision by majority of repeated measurements. If 
the state of the target qubits is the eigenstates of the c-U, the unitary transform 
results in the same phase value to the control qubit. However, in most case, the 
target qubit state is a superposition of the eigenstates. In this case, the total 
state at the k th operation step after the rotation gate is transformed by the 
Hadamard gate to 



1 



1 



1 



|0> 



1 



K), (16) 



where the states \uq) and jitf) defined by 

ko) = C'o ^ Cs \Us) 

Wl) =Ci Cs \Us) 



(17) 



are normalized as to {uq\u'^) + {u\\u\) = 1. Suppose measurement outcome is 
"0", then the state collapse to 



1 



:Mo|*) 



1 



where Pq = [1 + (2 



1 



V 



\^)\uo) 



1 - V 



1) \ui) 



(18) 



l)wcosi5] /2 is the probability to obtain the out- 



come "0", and the states |mo) and are defined by 



1 



1- e' 



"1 



K> + 



AS 



(19) 



Partial trace on the control qubit yields a density matrix of the target qubit 
states 



1 



, , ,-„.-„ m 

The c-U and Hadamard transform H16() followed by the measurement (|13|) on 
the density matrix (|20|) results in the same value of the error probability (|15|) . 
Therefore, the accumulation of M measurement outcomes will reduce the error 
probability to 

PM^ E [ )P''-Hl-P)\ (21) 

j=0 ^ 

which decreases rapidly for p <C 1. For example, the error probability per 
operation p — 0.07 and accumulation M = 10 yield the error probability pio = 
3 X 10^*, which agrees with the value obtained from the experiment on 1024 
qubits. 
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4.3 Toward quantum computers 

The main drawback of the MQFT is time, because it operates seriahy. The 
response time of the current circuit is hmited by electronics. We may expect 
the response time to be reduced to several nanoseconds. If the operation time 
is decreased to 1 ns, the target qubits should remain coherent for at least 10 jis 
to complete controlled-unitary operations with a thousand control qubits (1 ns 
X 1000 qubits x 10 accumulations.) This coherence time would be possible, be- 
cause the coherence time between the meta-stable states of an atom reaches 0.5 
ms^Sl- Here, we consider a quantum computer consists of photon control qubit 
and atom target qubits. The input photons induce a unitary transformation on 
the atom qubits. The controlled operation can be achieved by the polarization 
selection rules in the atomic transitions, for example. The photons are scattered 
by the atoms and obtain a phase shift through the kick-back effect. The phase 
shift will be analyzed by MQFT and provide a solution of the problem. This 
scheme shows an interesting analogy between the quantum computation and 
the spectroscopy; the change in the system (target qubits) state is probed by 
the change in the scattered photon state (control qubit.) 

The remaining problems are in realizing controlled-unitary operations; the 
QFT by itself will not provide an exponential speed up in comparison with 
classical algorithms. One possibility is using atomic target qubits as discussed 
above, where the unitary transformation is achieved by the interaction between 
atoms. We need to construct unitary transformations in the atomic qubits, so 
that it may be seen just moving the difficulty to the atoms. Nevertheless, we 
believe this scheme would make the problem simpler. It can exploit fairly large 
interaction between atomic qubits to obtain two qubit gates. Measurement on 
photonic control qubit will resolve the read-out problem. Another possibility 
is to combine with the linear-optics gates proposed by Knill, Laflamme, and 
MilburnfT^. Since the KLM scheme utilizes single photon states, it suits well 
to the present MQFT circuit. 

The present MQFT circuit may find its applications on such as finding eigen- 
values and eigenvectors [T5] and clock synchronization^^ in the near future. In 
these applications, fewer controlled-unitary operations are enough to achive a 
meaningful task than those required in the factorization algorithm. 

5 CONCLUSION 

We have shown successful MQFT using single photon by the use of commer- 
cially available fiber optic devices. The MQFT is robust to errors as suggested. 
We have also shown that decision by majority is simple to implement and is 
useful to reduce error probability. Our implementation can be used in the fi- 
nal part of the order-finding algorithm. In other words, the 'easier part' of 
the Factorization algorithm is ready. The present MQFT circuit decreases the 
qubit-gate-operation product, and thus reduces the requirement for the thresh- 
old on the error probability. We believe the present MQFT implementation 
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provides a useful base to realize a large scale quantum computation. 
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